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ABSTRACT 

The injection of ultra-high energy cosmic rays in the intergalactic medium leads to the production of a GeV-TeV gamma-ray halo 
centered on the source location, through the production of a high electromagnetic component in the interactions of the primary 
particles with the radiation backgrounds. This paper examines the prospects for the detectability of such gamma ray halos. We explore 
a broad range of astrophysical parameters, including the inhomogeneous distribution of magnetic fields in the large scale structure as 
well as various possible chemical compositions and injection spectra; and we consider the case of a source located outside clusters of 
galaxies. With respect to the gamma-ray flux associated to synchrotron radiation of ultra-high energy secondary pairs, we demonstrate 
that it does not depend strongly on these parameters and conclude that its magnitude ultimately depends on the energy injected in the 
primary cosmic rays. Bounding the cosmic ray luminosity with the contribution to the measured cosmic ray spectrum, we then find 
that the gamma-ray halo produced by equal luminosity sources is well below current or planned instrument sensitivities. Only rare 
and powerful steady sources, located at distances larger than several hundreds of Mpc and contributing to a fraction > 10 % of the flux 
at 10" eV might be detectable. We also discuss the gamma-ray halos that are produced by inverse Compton/pair production cascades 
seeded by ultra-high energy cosmic rays. This latter signal strongly depends on the configuration of the extragalactic magnetic fields; 
it is dominated by the synchrotron signal on a degree scale if the filling factor of magnetic fields with B > 10^'* G is smaller than a 
few percents. Finally, we discuss briefly the case of nearby potential sources such as Centaurus A. 
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1. Introduction 

The quest for the sources of ultrahigh energy cosmic rays has 
long been associated with the search of their secondary radiative 
signatures. While propagating, the former indeed produce very 
high energy photons through the interactions with the ambient 
backgrounds; and gamma rays should be valuable tools for the 
identification of the birthplace of the primary cosmic rays as they 
travel in a straight manner, contrarily to charged particles. 

The detection of such photon fluxes is far from straightfor- 
ward, however. On purely astrophysical grounds, the propaga- 
tion of gamma rays with energy exceeding several TeV is ob- 
structed by their relatively short pathlength of interaction with 
cosmic microwave background (CMB) and infrared photons. 
These interactions lead to the production of high energy electron 
and positron pairs which in turn up-scatter CMB or radio pho- 
tons by inverse Compton processes, initiating electromagnetic 
cascades. One does not expect to observe gamma rays of en- 
ergy above ~ 100 Te V from sources located beyond a horizon of 
a few megaparsecs ( Gould & Schredei 1967 : Wdowczvk et alj 



ll972HProtheroell986l:lProtheroe & Stanevll993h . Below fliisen- 
ergy, the detection of the gamma ray signal originating from ul- 
trahigh energy cosmic rays depends on the luminosity and angu 
lar extension of the observed object and obviously on the sensi- 
tivity and angular resolution of the available instruments. 
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Intergalactic magnetic fields inevitably come into play in this 
picture: i) primary cosmic rays can be deflected by the mag- 
netic field surrounding the source prior to the production of sec- 
ondary photons/electrons/positrons, ii) secondary electrons and 
positrons produced by cosmic ray interactions or by photon- 
photon interactions can also be deflected by the intergalactic 
magnetic field during the cascade. In both cases the ultimate 
gamma ray emission can experience a sizable spread around the 
ultrahigh energy cosmic ray source. 

lAharonianl (|2002|) and lGabici & Aharoniar] (|2005|) have ob- 
served that the magnetic field environment of the source may 
also lead to the generation of a multi-GeV gamma ray halo 
around ultra-high energy cosmic ray sources. The idea is the fol- 
lowing: for high enough energies and in the presence of a mag- 
netic field large enough in the surroundings of the source, the 
electron and positron pairs produced through the interaction of 
ultra-high energy particles with the ambient backgrounds may 
lose a significant fraction of their energy through synchrotron 
emission. This synchrotron component falls around ~ 10 GeV, 
below the pair creation threshold, hence it will not be affected 
by further electromagnetic cascading. Assuming that the source 
is located at 100 Mpc distance, at the center of a homoge- 
neously magnetized sphere of radius 20 Mpc with B = 1 nG, 
iGabici & Aharonian (2005) estimate that this si gnal would have 
an an gular size of a fraction of degree (see also Aharonian et al.l 
I2010h . It could therefore appear point-like for instruments such 
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as the Fermi Space Telescope but as an extended source for fu- 
ture imaging Cerenkov telescopes. 

This process is particularly interesting, because its signature 
could easily be differentiated from a point-like gamma ray signal 
emitted by leptonic or hadronic channels inside the source, asso- 
ciated with the acceleration of particles to energies possibly well 
below 10'"^ eV. The dominant background to this synchrotron 
halo is actually the halo produced by multi-TeV electrons that 
undergo inverse Compton on the cosmic microwave background 
and that were se eded by multi-TeV ph otons at > Mpc distance 
from the source dAharonian et al.l[l994l) . However, as discussed 
in the present paper, discrimination between these two signals 
should be possible thanks to the different dependencies of the 
angular images on physical parameters and on energy. Then, 
one could argue that the detection of a synchrotron halo around 
a powerful source would constitute an unambiguous signature 
of acceleration to ultra-high energies in this source: first of all, 
the inverse Compton cooling length of electrons with energy 
1 TeV < Ef < 10'^ eV is smaller than 300 kpc, hence the 
gamma-ray signal of very high energy electrons produced in the 
source would appear point-like at distances > 100 Mpc; further- 
more, secondary electrons generally carry a few percents of the 
parent proton energy, so that > 10'^ eV secondary electrons cor- 
respond to > 10'^ eV protons; finally, if the radiating electrons 
are seeded away from the source by high energy nuclei, these 
nuclei must carry an energy > lO'^eV as lower energy nuclei 
are essentially immune to radiative losses. 

Given our state of knowledge on the sources of ultrahigh en- 
ergy cosmic rays, the detection of such a halo would have a last- 
ing impact on this field of research. The detectability of such 
signals thus deserves close scrutiny. 

The production of gamma ray signatures of ul trahigh energy 
cosmi c ray s have been stud ied numerically by iFerrigno et alj 
(l2004l) and lArmengaud et~a l. (2006). However, the former au- 
thors have focused their study on the Compton cascades follow- 
ing the production of ultrahigh energy photons and pairs close to 
the source and they have neglected the deflection imparted by the 
surrounding magnetic fields which dilutes the gamma-ray sig- 
nal (Gabici & Aharonian 2005), while Armengaud et al. (2006) 
have addressed both the synchrotron emission of secondary pairs 
and the Compton cascading down to TeV energies, albeit for 
the particular case of a source located in a magnetized cluster 
of galaxies. The high magnetic field that prevails in such envi- 
ronments increases the residence time of primary and secondary 
charged particles and thus also increases the gamma ray flux. 

The present paper aims at examining the prospects for the de- 
tectability of gamma ray halos around ultrahigh energy cosmic 
ray sources, relaxing most of the assumptions made in the above 
previous studies. In particular, we discuss the more general case 
of a source located in the field, outside clusters of galaxies. We 
focus our discussion on the synchrotron signal emitted by sec- 
ondary pairs, which offers a possibility of unambiguous detec- 
tion; nevertheless, the deflection and the dilution of the Compton 
cascading gamma r ay signal a t TeV e nerg ies is also discussed. 
Going further than lAharonianl (l2002h and iGabici & AharonianI 
(l2005h . we take into account the inhomogeneous distribution of 
the magnetic fields in the source environment. We also relax the 
assumption of a pure proton composition of ultra high energy 
cosmic rays, underlying to the above studies. The chemical com- 
position of ultrahigh energy cosmic rays indeed remains an open 
question. While experiments such as the Fly's Eye and HiRes 
have suggested a transition from heavy to l ight above ~ lO'^'^ eV 
(iBird et alJl993l;lAbbasi et alJ2005Ll20Toh . the most recent mea- 
surements made with the Pierre Auger Observatory rather point 



towards a heavy cornposition above lO'^ e V (lUnger et al.ll2007l : 
[Abraham et al.l2009l : lAbraham et al.l2010l) . As the energy losses 
and magnetic deflection of high energy nuclei differ from those 
of a proton of a same energy, one should naturally expect differ- 
ent gamma ray signatures. 

The lay-out of the present paper is as follows. In Section|2] 
we first test the dependence of the gamma ray flux produced by 
ultrahigh energy cosmic rays on the type, intensity and struc- 
ture of magnetized environments. We also discuss the effects of 
various chemical compositions and injection spectra. We con- 
clude on the robustness of the gamma ray signature according 
to these parameters and find that the normalization and thus the 
detectability of this flux ultimately depends on the energy in- 
jected in the primary cosmic rays. In Section [3] we discuss the 
detectability of ultrahigh energy cosmic ray signatures in gamma 
rays. Applying the results of our calculations, we show that the 
average type of sources contributing to the ultrahigh energy cos- 
mic ray spectrum produces a gamma ray flux more than two or- 
ders of magnitudes lower than the sensitivity of the current and 
upcoming instruments. We then explore the case of rare pow- 
erful sources with cosmic ray luminosity over energy 10'"^ eV 
of Li9 > 10"*'*"'*^ erg s"'. We assume throughout this paper that 
sources emit isotropically and discuss how the conclusions are 
modified for beamed emission in SectionH] The gamma ray sig- 
natures of those sources could be detectable provided that they 
are located far enough not to overshoot the observed cosmic ray 
spectrum. Finally, we also briefly discuss the detection of nearby 
sources, considering the radiogalaxy Centaurus A as a protoypi- 
cal example. We draw our conclusions in Section|4] 

2. Simulations 

As mentioned above, we focus our discussion on the synchrotron 
signal that can be produced close to the source by very high en- 
ergy electrons and positrons, that result themselves from the in- 
teractions of primary cosmic rays. The signal associated to in- 
verse Compton cascades of these electrons/positrons on radia- 
tion backgrounds will be discussed in Section [33l 

The secondary electrons and positrons are created through 
one of the three following channels: i) by the decay of a charged 
pion produced during a photo-hadronic interaction {Ay ^ + 

then n'^ — > yu^ v^, and //^ — » v^, with A a cosmic ray 
nucleus), ii) by photo pair production during an interaction with 
a background photon {Ay — » + ...), or iii) by the disin- 

tegration of a neutral pion into ultrahigh energy photons which 
then interact with CMB and radio backgrounds to produce elec- 
tron and positron pairs {Ay — » + then if — > 2y, and 
yybg e^e", with ybg a cosmic background photon). In all 
cases, the resulting electrons and positrons typically carry up to 
a few percents of the initial cosmic ray energy. 

While propagating in the intergalactic medium, these ul- 
trahigh energy electrons and positrons up-scatter CMB or 
radio photons through inverse Compton processes and/or 
they lose energy through synchrotron radiation. Following 
iGabici & AharonianI ( l2005h . the effective inverse Compton cool- 
ing lengtlfl on the CMB and radio backgrounds can be written 



' Let us recall that at very high energy, pair production (yybg — > 
e+e") transfers energy to one of the pairs. The inverse Compton scatter- 
ing of photons (eyhs — > ey) occurring in the Klein-Nishina regime, 
nearly all the electron or positron energy is again transferee! to the 
up-scattered photon. Thus the initial electron energy in a cascade is 
degraded very slowly, until the electron or positron energy is either 
radiated in synchrotron, or the photon energy falls beneath the pair 
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as x,y K SMpcC^JlO'^^eVy'S with a^c = 1 if the electron 
energy E, < 10'** eV and q-k ^ 0.25 if 10"^ eV < < 10^" eV. 
Above 10'** eV, the scaling of the cooling length actually de- 
pends on the assumptions made for the radio background, which 
is unfortunately not very well known. The possible differences 
that such uncertainties could introduce should however not af- 
fect the results discussed in this paper. This length scale has to 
be compared with the synchrotron cooling length 



'■'''P'^(io^) (lo^) 



(1) 



where B is the magnetic field intensity (assumed homogeneous 
over this distance) a nd Eg the electron or positron energy. 

As discussed in iGabici & AharonianI (l2005h . the opposite 
scalings of x^y and x^b with electron energy imply the exis- 
tence of a cross-over energy E^, which depends on B and is such 
that beyond E^, electrons mainly cool via synchrotron instead 
of undergoing an inverse Compton cascade. For a magnetic field 
of intensity B = 10 ''fincG, fix ~ lO^^eVfi-^!. for B„g > 1, 



10"* eVB:A'^ for 1 > finC > 0-1 and £x 



.^. .^ ^nc, ^v... ....^^x - 10^"eVfi^-' for 

0.1 > BgG, see Fig. 1 o f iGabici & AharonianI (l2005h and seelLed 
(Il998h for further details on the cascade physics. 

Finally, the emitted synchrotron photon spectrum peaks at 
the energy: 



£y,sy„ ~ 68 GeV 



llOnG/UOi^eV/ 



(2) 



2.1. Magnetic field configuration 



Our current knowledge on the structure of magnetic fields at 
large scales is very poor. It mainly stems from the lack of ob- 
servations due to the intrinsic weakness of the fields. Moreover, 
no satisfactory theory has been established to explain the ori- 
gins of these fields and the further magnetic enrichment of the 
Universe. There is now a decade of history of numerical simula- 
tions that endeavor to model the distribution of the ma gnetic field 
at larg e scales. One may cit e in p articular the works of iRvu et"aD 
(ll998i) JDolag et al. (2004, 2"005l) . ISigl et"al] ( |2004 and most re- 
cently |Das_etalj (12008). In these studies, magnetized seeds are 
injected in a cosmological simulation and the overall magnetic 
field then evolves coupled with the underlying baryonic matter 
The final amplitude of the field is set to adjust the me asured in- 
tensity inside clusters of galaxies Cm iDas efal 1 120081 the ampli- 
tude is not arbitrary but is directly estimated from the gas kinetic 
property and corresponds nevertheless to the observations inside 
clusters). Unfortunately, these various simulations do not con- 
verge as to the present-day configur ation of extragalactic mag- 
netic fields (see for example Fig. 1 of lKotera & Lemoinell2008al) 
and consequently, they also diverge as to the effect of these mag- 
netic fields on the transport of ultrahigh energy cosmic rays. The 
cause of this difference likely lies in the choice of initial data, 
which cannot be constrained with our current knowledge. Given 
the fact that such simulations are time consuming, these tools 
cannot be considered as practical. 

In view of this situation, we explore the influence of three 
typical types of magnetic field structures on t he gamma ray emis- 
sion, following the method developed by iKotera & Lemoing 



production threshold. For this reason one may consider the cascade 
(e — > 7 — » e...) as the disintegration o f one single particle, that loses en- 
ergy over an effective loss length x^^ ( ISteckerlfl973l : lGould & Rephaelil 
[1971 . 



(I2008ah . We map the magnetic field intensity distribution ac- 
cording to the underlying matter density p, assuming the rela- 
tion B - Bi) fip), where Bq is a normalization factor, and the 
dimensionless function / may be modelled as: 



/iso(p) 
/ani(p) 

/contrast(p) — P 



^2/3 



^0.9 



1+1^ 

<P> 



-2 



n-1 



(3) 
(4) 

(5) 



Note that Bo corresponds approximately to the mean value of 
the magnetic field in the Universe. These phenomenological re- 
lationships greatly simplify the numerical procedure to obtain a 
map of the extragalactic magnetic field since a pure dark mat- 
ter (i.e. non hydrodynamical) simulation of large scale struc- 
ture provides a sufficiently good description of the density field. 
These relationships are motivated by physical arguments related 
to the mechanism of amplification of magnetic fields during 
structure formation and they encapture the scalings observed 
in numerical simulations up to the uncertainty inh erent in such 
sim ulations (see the correspo nding discussion in iDolad l2006l 
and lKotera & Lemoinel2008ah . The first relation is typically ex- 
pected if baryonic matter undergoes an isotropic collapse. In 
dense regions of the Universe where viscosity and shear effects 
govern the evolution of the magnetized plasma however, the rela- 
tion would get closer to Eq. (|4]i. This formula can be derived an- 
alytically in the case of an anisotropic collapse al ong one or two 
dimensions, in filaments or sheets for example dKing & ColesI 
I2OO6I) . The last model is an ad-hoc modeling of the suppression 
of magnetic fields in the voids of large structure which leaves 
unchanged the distribution in the dense intergalactic medium 
(meaning p > {p}). It thus allows to model a situation in which 
the magnetic enrichment of the intergalactic medium is related 
to structure formation, for instance through the pollution by 
starburst gala xies and/or radio-galaxies (see the corresponding 
discussion in iKotera & Lemoinell2008bl) . In the following, we 
will refer to these three models as "isotropic", "anisotropic" and 
"contrasted" respectively. 

In order to have a representative sample of the underlying 
matter density at large scales, we use a three-dimensional out- 
put (at redshift z = 0) of a cosmological Dark Matter simula- 
tion provided by S. Colombi. It assumes a ACDM model with 

= 0.3, Q.A = 0.7 and Hubble constant h = //o/(100 km s"' 
Mpc"') = 0.7. The simulation models a 200/1"' Mpc comoving 
periodic cube split in 256^ cells, where the Dark Matter overden- 
sity is computed. We do not resolve structures below the Jeans 
length, which implies that we can identify the computed Dark 
Matter distribution to a gas distribution. 

2.2. Numerical code 

We compute the gamma ray signal through num erical Monte 
Carlo simulations, using the code presented in IKotera et al.l 
(I2OO9I) . In a first step, we propagate ultrahigh energy cosmic 
rays (protons and heavier nuclei) in one of the three-dimensional 
magnetic field models described above. The coherence length 
Ab of the magnetic field is assumed to be constant throughout 
space and is set to 100 kpc. This certainly is somewhat ad- 
hoc but: the coherence scale enters the calculation in the de- 
flection angle under the form B yfAs, hence discussing various 
normalizations of the magnetic field strength, as we do further 
below, allows to encapture different possible values for Ab\ fur- 
thermore, 100 kpc seems a reasonable compromise between an 
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Fig. 1. Distribution of the magnetic field intensity in the region surrounding an ultrahigh energy cosmic ray source. Each panel is a 
1 . 1 Mpc thick (size of a grid cell) slice of the simulated Universe, cut along the axis perpendicular to the observed plan. From left to 
right, the distribution is of type "isotropic" (Eq.O, "anisotropic" (Eq.|4| and "contrasted" (Eq.|5]l. The color code is the same for all 
the panels and the values indicated on the right hand side correspond to the logarithm of the magnetic field intensity in nano Gauss. 
The black line indicates the axis of the filament that is studied throughout this paper and the black square the chosen position of the 
source. The normalization factor equals Bq - I nG. 



upper limit of a few hundreds o f kpc for Ab based on the turn- 
around time of largest eddies (IWaxman & Bahcalil 1999h and 
the typi cal galactic scales ~ 10- 30 kpc, see Kotera & Lemoin j 
(l2008ah for a detailed discussion. The generation of secondary 
photons and electrons and positron s through pion prod uction is 
treated in a discrete manner, as in lAllard et all ( |2006|) for nu- 
clei p rojectiles with A > 1 and with SOPHIA dMiicke et"an 
1 19991) for proton and neutron cosmic rays. Electron and positron 
pair c reation through photo-ha dronic processes is implemented 
as in lArmengaud et al.l (l2006h . For each time step (chosen to 
be much larger than the typical mean free path of pair photo- 
production processes), we assume that an ensemble of electron 
and positron pairs are generated with energies distributed ac- 
cording to a power law, and maximum energy depending on the 
primary particle. The use of an impro ved pair spectrum as com- 
puted in iKelner & AharonianI (l2008h would lower the gamma 
ray flux resulting from direct pair production by a factor of a 
few in the range Ey < lOGeV, and leave unchanged the pre- 
diction at a higher energy (Armengaud, private communication). 
Furthermore, we will show in the following that direct pair pro- 
duction processes provide a sub-dominant contribution to the 
overall gamma ray flux in the range Ey > 0.1 GeV, hence the 
overall diff'erence is expected to be even smaller. 

We investigate in this paper the possible detection of photons 
in the GeV-TeV energy range. Equation ^ suggests that, unless 
the source is embedded in a particularly strong field, the major 
contribution in this range will come from electrons of energy 
Ef > 10"^ eV. The flux will thus essentially result from cosmic 
rays of energy higher than E ~ 10'^ eV. For this reason, we 
chose to inject cosmic rays at the source between E,^in - 10'^ eV 
et -Emax = 10^"'^ eV. We assume our source to be stationary with 
(isotropic) cosmic ray luminosity integrated over E = 10'*^ eV 



of L 



'£,19 = 10"*^ erg s ', and a spectral index of = 2.3. 



We take into account photo-hadronic interactions with CMB 
and infrared photons. The diffu se infrared ba ckground is mod- 
elled according to the studies of lStecker et alj (|2006^. We do not 
take into account redshift evolution, as its effect is negligible as 
compared to the uncertainties on our other parameters. We also 
neglect baryonic interactions in view of the negligible density in 
the source environment. 



Once the secondary particles produced during the propaga- 
tion have been computed, we calculate the synchrotron pho- 
ton fluxes, taking into account the competition with inverse 
Compton scattering by Monte Carlo over steps of size 100 kpc. 
The mean effective electro n loss length is calculated following 
Gabici & AharonianI ( |2005|) . using the radio background data of 
Clark etalJ (119701) . We assume that each electron emits a syn- 
chrotron photon spectrum dNy/dEy oc Ey^^^ with a sharp cut-off 
above ~ (2/3) -Ey.syn, where fiy.syn is given by Eq. 

Regarding the photons produced through the neutral pion 
channel, we draw the interaction positions with C MB and ra - 
dio photons using the mean free paths computed bv lLeel(ll998h . 
We assume that one of the electron-positron pair inherits of the 
total energy of the parent photon. 

2.3. Synchrotron signal for various magnetic field 
configurations 

We first study the influence of inhomogeneous magnetic fields 
on the gamma ray flux produced by ultrahigh energy cosmic 
rays near their source. In particular, we examine the case of a 
source placed in a rather dense region of a filament of large 
scale structure, see Fig.[T] On the transverse scale of the structure 
l^, of order of a few Mpc, one expects the secondary electrons 
and positrons to be distributed isotropically around the source. 
A population of ultrahigh energy cosmic rays seeds the inter- 
galactic medium with secondary pairs in a roughly homogeneous 
manner on the energy loss length scale of the interaction process 
A, i.e. the number of electrons/positrons injected per unit time 
per unit distance to the source dNe/drdt does not depend on r 
as long as r «c /I. For ultrahigh energy protons, this length scale 
is ~ 1 Gpc above 10''^eV for pair production and ~ 100 Mpc 
above 10"" eV for pion production. From the point of view of 
secondary synchrotron y emission, however, the emission region 
is limited to the fraction of space in which the magnetic field is 
sufficiently intense for synchrotron cooling to predominate over 
inverse Compton cooling. 

The bulk of gamma rays is produced by electrons of energy 
Ee ^ 10'^ eV, themselves produced through pion production 
by protons of Ep > 10^'^ eV. In the absence of significant de- 
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Fig. 2. Flux of photons produced by synchrotron emission for a filament observed from the directions forming the indicated angles 
with its axis. The filament is located at 100 Mpc and harbors a source that injects a pure proton composition with (isotropic) 
luminosity of 19 - 10"^^ erg/s and a spectral index of 2.3. The four panels represent various models of magnetic fields. From 
left to right: "anisotropic" with Bo = 0.1 nG, "anisotropic" with Z?o = 1 nG, "isotropic" with Z?o = 1 nG and "contrasted" with 
Bo = 1 nG. 



flection experienced by these primary particles, one expects the 
synchrotron emission to be enhanced along the filament axis, 
where the magnetic field keeps a high value over a longer length 
than along the perpendicular direction (see Fig. [TJ. The whole 
structure where x^b < Xey, i-e. -Bigm <: 0.5nG(£'j,/10'^eV)"-'^'' 
is then illuminated in synchrotron from electrons and positrons 
provided x^b < cfesc, where fgsc denotes the escape timescale 
from the magnetized region. For a magnetic field with B ~ 
lOnG and coherence length Ab < 100 kpc, the Larmor time 
fL ~ lMpc(£',/10'''eV)(B/10nGr' ~ ^b, hence confine- 
ment is marginal and the escape length fgsc ~ l± '» x^b (see 
ICasse et alj|2002l for more details on the escape length). 

At a given energy E^, the ratio between the fluxes along and 
perpendicular to the filament should be of order a/b, where a 
and b are the characteristic lengths of the axis of the filament 
where Bigm > 0.5 nG£';^^'*. 

Figure |2] depicts indeed such an effect. It represents the 
synchrotron emission produced by a filament at a distance of 
100 Mpc, embedding a source of luminosity L^ ig = lO''^ erg s"' 
and an injection spectral index of 2.3. One can notice that the dif- 
ference between the highest and lowest flux values is of a factor 
~ 2 - 5, which corresponds to the ratio between the filament ma- 
jor axis. For a narrower or longer filament, this interval would 
be more pronounced - which does not necessarily mean that the 
flux would be enhanced along the filament axis. 

This diff'erence of flux according to the observation angle 
disappears at low energy for the second and third panels of 
Fig. |2] This is intricately related to the combined effect of the 
configuration of the magnetic field and the ratio of electrons pro- 
ducing the observed peak of gamma rays. It corresponds to a co- 
incidence and should not be viewed as an effect of diffusion of 
primary protons. We tested indeed that these differences remain 
present when cosmic rays are forced to propagate rectilinearly. 
For the contrasted magnetic field, the gamma-ray fluxes are over- 
all lower, as the field strength diminishes sharply outside of the 
central part of the filament. The orthogonal section of the fila- 
ment is thus substantially magnetized only over a short length, 
which consequently reduces the gamma-ray flux observed from 
90°. For the isotropic and anisotropic cases with Bo = 1 nG, elec- 
trons and positrons experience a weaker field on average when 
they cross the filament orthogonally, instead of traveling through 



their entire length. For the flux observed at 90°, photons with en- 
ergy < 10 GeV are thus generated by electrons of slightly higher 
energy (£90) than in the 0° case (Eq). By coincidence, £90 cor- 
responds approximately to the peak of distribution of secondary 
electrons, whereas E^ is shifted according to the peak. The flux 
observed at 90° is thus amplified for energies < 10 GeV due to 
the enhanced number of electrons contributing to the emission, 
which brings up the fluxes at various observation angles closer. 
Finally, for the anisotropic case with Bo = 0.1 nG, the emission 
is suppressed when the filament is not observed along its major 
axis, because of the weak mean field that particles experience 
(the mean field remains reasonably strong along the axis). The 
gap between the flux at 0° and the other cases stems from this 
effect. 

Most of all, Fig.|2]demonstrates that the photon flux around 
Ey ~ 10 GeV is fairly robust with respect to the various mag- 
netic field configurations, which can be understood from Fig.[T] 
the extent of the magnetized region above B ~ 1 nG is indeed 
not fundamentally modified by the chosen distribution and the 
photon flux is hardly more sensitive to the overall intensity of 
the field (in particular, we tested that a value of Bo = 10 nG does 
not affect the results of more than a factor 2, for all magnetic 
configurations). All in all, for a quite reasonable normalization 
of the magnetic field, meaning an average Bo ~ 0.1 - lOnG at 
density contrast unity, and for various scalings of the magnetic 
field with density profile, the variations in flux remain smaller 
than an order of magnitude in the range Ey ~ 10 - 1000 GeV. 

Figure |2] also illustrates that the energy at which the signal 
peaks depends only weakly on the normalization of the mag- 
netic field Bo, and this remains true up to a few nano-Gauss. This 
can be understood qualitatively, assuming for simplicity the field 
to be homogeneous over the area of interest, with strength Bo. 
For Bo = 1 nG, the electron cross-over energy Ex beyond which 
electrons radiate in synchrotron reads E^ - lO'^ eV (see the dis- 
cussion in Sec. |2]i. However, the energy spectrum EgdNe/dEe 
of secondary electrons deposited in the source vicinity peaks 
at ~ 10'^ eV, as a result of the competition between the oppo- 
site scalings with energy of the abundance of parent protons 
and of the parent proton energy loss rate. Therefore, the syn- 
chrotron signal is expected to peak at Ey ~ 10 GeV, see Eq. |2] 
with Bo = 1 nG and E, = 10''' eV. Now, for Bo = 0.1 nG, the 
cross-over energy Ex - 4- x 10'^ eV so that one should use 
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E, = 4x 10'''eV in Eq.|2]with Bq = 0.1 nG, giving Ey ~ lOGeV 
again: the larger typical electron energy (among those radiating 
in synchrotron) compensates for the smaller value of Bq. As Bq 
is increased significantly above 1 nG, the peak location of the 
gamma-ray signal will increase in proportion, as most of the sec- 
ondary electron flux can be radiated in synchrotron. 

Note that the flux should be cut off above Ey ~ 10 TeV due 
to the opacity of the Universe to photons above this energy. This 
effect is not represented in these plots for simplicity, as its pre- 
cise spectral shape depends on the distance d to the source (while 
the sub-TeV gamma ray flux scales in proportion to L^Jcf-). 

At lower energy {Ey < 1 GeV), the flux intensity can vary 
by more than two orders of magnitude depending on the chosen 
magnetic configuration. In the low energy range, magnetic con- 
finement indeed starts to play an important role. Furthermore, 
the inverse Compton energy loss length shortens drastically at 
these energies so that only strong magnetic fields can help avoid- 
ing the formation of an electromagnetic cascade and lead instead 
to synchrotron emission in the filament. 

2.4. Synchrotron signal for various cliemical compositions 
and injection spectra 

In this section, we discuss the effect of the injected chemical 
composition and of the spectral index at the source. As discussed 
above, there are conflicting claims as to the measured chemical 
composition at ultrahigh energies; in particular, the Hi Res exper- 
iment reports a light composition dAbbasi et al while the 
Pierre Auger Observatory measurements point toward a compo- 
sition that becomes increas ingly heavier at energies above the 
ankle ([Abraham et al.ll20Toh . Theory is here of little help, as the 
source of ultrahigh energy cosmic rays is unknown; protons are 
usually considered as prime candidates because of their large 
cosmic abundance, but at the same time one may argue that a 
large atomic number facilitates acceleration to high energy. 

In this framework, iLemoine & WaxrnanI (l2009l) have pro- 
posed a test of the chemical composition of ultra-high energy 
cosmic rays on the sky, using the anisotropy patterns measured 
at various energies, instead of relying on measurements of the 
depth of maximum shower development. It is shown in particu- 
lar that, at equal magnetic rigidities E/Z, one should observe a 
comparable or stronger anisotropy signal from the proton com- 
ponent than from the heavy nuclei component emitted by the 
source, even if the source injects protons and heavy nuclei in 
equal numbers at a given energy. When compared to the 99% 
c.l. anisotropy signal reported by the Pierre Auger Observator y 
at energies above 5.7 x 10''' eV ([Abraham et al.ll2007l l2008al) . 
one concludes that: if the composition is heavy at these energies, 
and if the source injects protons in at least equal number (at a 
given energy) as heavy nuclei, then one should observe a com- 
parable or stronger anisotropy pattern above 5.7 x lO'^/ZeV. 

Future data will hopefully cast light on this issue, but in 
the meantime it is necessary to consider a large set of possible 
chemical compositions. We thus consider the following injection 
spectra: 

1. Two pure proton compositions with indices a - 2.3 and 
2.7. A soft spectral index of ~ 2.7 seems to be favored to 
fit the observed cosmic ray spectrum in the case of a pure 
proton composition, especially in the "dip-model" proposed 
by Berezinsky et al. [200^, where the transition between the 
Galactic and extragalactic components happens at relatively 
low energy (E ~ 10'^'^ eV). 



2. A proton dominated mix composition with spectral index 
a = 2.3, based on Galactic cosmic ray abundances as in 
lAllard et a l. (2006). Such a spectrum is mostly proton dom- 
inated at ultra-high energies. 

3. A pure iron composition with spectral index a = 23. 

4. A mix ed composition that was proposed by lAllard et"al1 
(120081) . that contains 30% of iron. In this injection, the max- 
imum proton energy is Etaax.p - lO'** eV (jAllard et al. 200*81 
considered iSmax,/; = 4 x lO'^ eV but an equally satisfying fit 
ot the spectrum can be obtained for a lower value). Assuming 
that heavy elements are accelerated at an energy Z x fimax.p, 
and because of the propagation effects that disintegrate inter- 
mediate nuclei preferentially as compared to heavier nuclei, 
one would then measure on the Earth a heavy composition at 
the highest energy. The resulting detected spectrum allows 
to reproduce the maximum depth of shower measurements 
of the Pierre Auger Observatory. However, the anisotropy 
pattern should also appear in this case at energies of order a 
few EeV, as discussed above. An injection spectrum of 2.0 is 
necessary to fit the observed cosmic ray spectrum. 

For the first three cases, the maximum proton injection energy 
is ^max.p = 10^°-^ eV. In all cases, we assume E„,nx.z = Z x 
Em-dx.p for a nucleus of charge number Z (an exponential cut-off 
is assumed). 

Figure |3] presents the fluxes obtained for these different in- 
jection spectra. The source is again located in the filament de- 
scribed earlier, at a distance of 100 Mpc, and has a luminosity 
above E = lO'** eV of L^jg = 10''^ erg s"'. Due to the compar- 
atively larger number of protons at low energy, the photon flux 
below 1 GeV is amplified in the case of an injection spectrum of 
2.7 as compared to a 2.3 index. For similar reasons, the flux at 
higher energies (> 10 GeV) is lower. 

The shift in amplitude between the fluxes produced by the 
different compositions can be understood as follows. The energy 
of the electrons and positrons produced by photo-disintegration 
is proportional to the nucleus Lorentz factor £(. oc oc Ea/A. 
Hence at a given electron energy, the amount of injected en- 
ergy by heavy nuclei is reduced if the spectral index is > 2. 
Besides, pion production via baryonic resonance processes on 
the CMB photons happens at energies higher than for pro- 
tons (E > 10^" eV). Obviously, the flux produced by the mix 
"Galactic" composition with a 2.3 injection differs only slightly 
from that produced by protons with the same spectral index as 
such a composition is dominated by protons. 

The above trends can also be noted in Fig. |4] which re- 
veals a smaller contribution of photo-disintegration processes 
(in black dotted and pale blue dash-dotted lines) than for a pro- 
ton composition. The yield of the electrons and positrons pro- 
duced through photo-pair production (in blue dashed lines) is 
not as strongly attenuated however, as the secondary nucleon are 
emitted above the pair production threshold (but below the pion 
photo-production threshold). The presence of a second bump at 
lower energy in the neutral pion channel (pale blue dash-dotted) 
is related to the production of photons via the giant dipole reso- 
nance process. 

Finally, the flux is suppressed in the case of an iron enriched 
mix composition, because of the relatively low cut-off energy, 
Emax.p - 10''' eV. In this case, /imax.p is lower than the pion pro- 
duction threshold energy for protons and the heavy nuclei maxi- 
mum energy is also lower than the pion production threshold, so 
that this channel of secondary e^e^ injection disappears. 
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Fig. 3. Photon flux produced by synchrotron emission for dif- 
ferent source injections. We chose an "anisotropic" distribution 
of the magnetic field with normalization Bq = 1 nG. The source 
has a luminosity of L^jg = 10"*^ erg s"', and is located in the 
filament described previously, at a distance of 100 Mpc. The av- 
erage flux integrated over aU angular lines of sight are presented. 
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Fig. 4. Same as Fig.|3] but for the particular case of a pure iron in- 
jection with spectral index 2.3. In red solid lines, the total flux, in 
blue dashed lines the contribution of the photo-pair production, 
in black dotted lines the flux due to charged pion production via 
photo-disintegration processes and to neutron disintegration, in 
pale blue dash-dotted lines, the contribution of the neutral pion 
channel. 

3. Discussion on detectability 

3. 1. Synchrotron signal from equal luminosity sources 

Our study indicates that around ~ 10 - 100 GeV, the gamma 
ray flux produced by the propagation of ultrahigh energy cosmic 
rays is fairly robust to changes in composition and in the con- 



figuration of magnetic fields. The normalization of the flux ulti- 
mately depends on the injected energy, namely, the cosmic ray 
luminosity and the injection spectral index at the source. The 
cosmic ray luminosity of the source may be inferred through 
the normalization to the measured flux of ultrahigh cosmic 
rays once the source density is specified. The apparent density 
of ultrahigh energy cosmic ray sources, which coincides with 
the true source density «s for steady isotropic emitting sources 
(which we assume here; see Section 4 for a brief discussion of 
bursting/beamed sources), is bound by the l ack of repeaters in 
the P ierre Auger dataset; «s > 10"^ Mpc""* jKashti & WaxmanI 
I2OO8I) . This lower bound may in turn be translated into an up- 
per limit for the gamma ray flux expected from equal luminosity 
cosmic ray sources, Lgjg < 10"*^ erg/s. 

Figure [5] presents the spectra obtained for sources of lu- 
minosity Le,i9 - 10"*" ergs"' and of number density = 
10"^ Mpc"-'. For each case, we represent the median spectrum 
of 1000 realizations of local source distributions, the minimum 
distance to a source being of 4 Mpc. Our calculated values are 
comparable to the cosmic ray fluxes detected by the two most re- 
cent experiments: HiRes and the Pierre Auger Observatory. The 
gamma ray fluxes obtained for single sources with this normal- 
ization (i.e. a luminosity of LE,\g - 10'^^ erg s"') at a distance of 
100 Mpc are presented in Figs. |2]|3] and |4] 
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Fig. 5. Propagated spectra of ultrahigh energy cosmic rays for 
various chemical compositions as described in section 12.41 and 
corresponding spectral indices to best fit the observational data. 
Sources are assumed to share an equal luminosity Le.i9 = 
10'*^ erg s"' and a spatial density = 10"^ Mpc"-'. 



These gamma ray fluxes lie below the limits of detectability 
of current and upcoming instruments such as Fermi, HESS or 
the Cerenkov Telescope Array (CTA), by more than two orders 
of magnitude. The sensitivity of the Fermi telescope is indeed of 



order ~ 2 x 10 GeV s cm around 10 GeV for a year of 
observation and the future Cherenkov Telescope Array (CTA) is 
expected to have a sensitivity of order ~ 10"" Ge V s"'cm"^ for 
100 hours for a point sour ce around TeV energies dAtwood et al.l 
l2009HWagner et al.ll2009h . 
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3.2. Rare and powerful sources 

We conclude from the previous subsection that only powerful 
sources embedded in magnetized environments, rare with re- 
spect to the average population of ultrahigh energy cosmic ray 
sources, may produce a gamma ray signature strong enough to 
be observed by current and upcoming experiments. In the fol- 
lowing, we examine the case of two such sources embedded 
in filaments, with (isotropic) luminosity chosen to be higher 
than the average allowed by the cosmic ray normalization for 
equal luminosity sources: Le,\9 - 10'*'* erg s"' and LEjg - 
10"*^ erg s"'. As a case study, we place these sources at a dis- 
tance of D = 100 Mpc and D - I Gpc respectively. Those dis- 
tances are consistent with the number density of objec ts found 
at com parable photon luminosities (see for example iWall et al.l 
I2005h . 



10^ 



10^ 



E 

> 



L,9=10''^ erg s"', D = 1 Gpc 

L,g=10" erg s"', D = 100 Mpc 

A Auger (PRL 2008) 



A A ^ 




19.0 19.5 
log E [eV] 



20.0 



20.5 



Fig. 6. Propagated spectra of ultrahigh energy cosmic rays from 
a single source with luminosity Lea9 - 10^^ erg s"' located at a 
distance of D = 1 Gpc (black solid) or L^ ig = 10'*'* erg s"' and 
D - 100 Mpc (green dashed). The injection s pectral index is 
2.0 in b oth cases. The observed Auger spectrum ([Abraham et al.l 
I2008bl) is overlaid. 



The ultrahigh energy cosmic ray spectra obtained after prop- 
agation from these two model sources are represented in Fig.|6] 
The injection at the source is assumed to be pure proton and 
with an index of 2.0. The comparison with the Auger data 
shows that the close-by source (green dashed) is marginally ex- 
cluded. More precisely, such a source should produce a strong 
anisotropy signal, of order unity relatively to the all-sky back- 
ground, if the magnetic deflecti on at this energy ~ 7 x 10 '^ eV 
is smaller than unity. Following lKotera & Lemoind (l2008bh . the 
expected deflection is of order 2°Z if one models the magnetized 
Universe as a collection of magnetized filaments with B = 10 nG 
{Ab = lOOkpc) immersed in unmagnetized voids with typical 
separation ~ 40 Mpc. This estimate neglects the deflection due 
to the Galactic magnetic field, which depends on incoming di- 
rection and charge. Qualitatively speaking, however, one would 
then expect the above source to produce a strong anisotropy sig- 
nal unless it injects only iron group nuclei with Z ~ 26 at the 
highest energies. In this latter case, however, the gamma-ray sig- 



nal would be lower by a factor ~ 5 or even suppressed if iSmax.p 
is less than the pion production threshold, as discussed in the 
previous Section. Furthermore, such a source might also give 
rise to a strong anisotropy at energies ~ 3 EeV, depending on 
the composition ratio of in jected protons to iron group nuclei 
dLemoine & Waxmanll2009l) . 

On the contrary, the remote source (indicated by the black 
solid line in Fig.|6l) contributes to ~ 10% of the total flux around 
Ed ~ 10*^ eV. Such a source would be wholly invisible in ultra- 
high energy cosmic rays, because of the large expected magnetic 
deflection at this energy, on a large path length of 1 Gpc. 

Let us now discuss the possibility of detecting these sources 
in gamma rays. Figure Q presents our computed images of syn- 
chrotron gamma ray emission for these two examples of sources 
embedded in the filament shown in Fig.[T] 

However, one should be aware that the sensitivity of gamma 
ray telescopes weakens for extended sources by the ratio of the 
radius of the emission to the angular resolution. Typically, Fermi 
LAT and HESS have an angular resolution of some fractions of 
a degree around energies of 10 GeV and 100 GeV respectively. 
CTA will have an even higher resolution of the order of the arc- 
minute above 100 GeV. For this reason, one will have to find a 
source with an acceptable range in luminosity and distance, so 
as to not lose too much flux due to the angular extension, but still 
being able to resolve this extension. Let us note that the magnetic 
configuration around the source should have a negligible role on 
the flux intensity, as we demonstrated in section |2] However, the 
magnetic configuration, in particular the strength and the coher- 
ence length, directly controls the extension of the image, which 
is given by the tran sverse displacement of the protons through 
magnetic deflection dGabici & Aharonianll2005h . 

Figure[8]demonstrates that the synchrotron gamma ray emis- 
sion from both our model sources could be observed by Fermi 
LAT and by CTA. Their flux integrated over the angular reso- 
lution of these instruments at ~ 10 GeV are indeed above their 
corresponding sensitivities. The case presented in dashed green 
line is however marginally excluded by the observed spectrum 
of cosmic rays and even more so by the search for anisotropy, as 
discussed previously (Fig.|6l). It thus appears that the gamma ray 
flux of a source of luminosity Le,i9 ~ 10'^* erg s"' will be hardly 
observable, as placing the source at a further distance to recon- 
cile its spectrum with the observations will dilute the gamma 
ray emission consequently, below the current instrument sen- 
sitivities. More promising are extremely powerful sources with 
Le,i9 = lO'*^ erg s"' located around 1 Gpc. According to Fig. [8] 
the emission from such sources would spread over a fraction of 
a degree, hence their image could possibly be resolved by Fermi 
and certainly by CTA. 

As discussed in the introduction, the observation of a multi- 
GeV extended emission as obtained in these images would pro- 
vide clear evidence for acceleration to ultra-high energies, pro- 
vided the possible background associated with a halo seeded by 



TeV p hotons can be removed. As discussed in [Aharoni an et al 
1 9941) . iDaietalJ (|2002|). Id'Avezac et alj (l2007i) . ,Elviv et al 



J2OO9I) and iN eronov & Semikod (l2009t) . the angular size of the 



latter signal is a rather strong function of the average inter- 
galactic magnetic field. In particular, if the field is assumed ho- 
mogeneous and of strength larger than 10"'"* G, the image is di- 
luted away to large angular scales. Accounting for the inhomo- 
geneity of the magnetic field (see Sec. l3.3l below for a detailed 
discussion), one concludes that the flux associated to this im- 
age, on a degree scale, is to be reduced by the filling factor of 
voids in which B < 10 '"^G within ~ 40 Mpc of the source 
(this distance corresponding to the energy loss distance of the 
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Fig. 7. Image of a filament in gamma rays produced via synchrotron emission, integrated over Ey - 1-100 GeV. We chose an 
"anisotropic" distribution of the magnetic field with normalization Bq - \ nG. Left: the filament is observed along its axis from a 
distance of 100 Mpc and the embedded source has a cosmic ray luminosity of L^jg = 10'*"^ erg s"'. Right: the filament is observed 
along its axis from a distance of 1 Gpc and the embedded source has a luminosity of L^jg - 10^^ erg s"'. The color bar indicates 
the logarithm of the flux in GeV cm - s White contours indicate a flux level of 10 '° GeV cm"^ s 



primary ~ 20TeV photons). Furthermore, the angular size of 
the halo seeded by multi-TeV photons depends strongly on pho- 
ton energy, i.e. 60 oc o r 59 oc /s^'^^ depending on Ab, see 
iNeronov & Semiko3 (120091) . In contrast, the synchrotron halo is 
more strongly peaked in energy around 10 GeV and in this range 
its angular size is essentially independent of energy, but is rather 
determined by the magnetic field strength and source distance. 
Assuming comparable total power in both signals, discrimina- 
tion should thus be feasible; of course discrimination would 
be all the more easier if the halo signal seeded by multi-TeV 
gamma-rays were diluted to large angular scales. 

The previous discussion thus indicates that detection might 
be possible for powerful rare sources at distances in excess of 
100 Mpc. However, one could not then argue, on the basis of 
the detection of a gamma ray counterpart that the source of ul- 
trahigh energy cosmic rays has been identified, as this source 
of gamma-rays must remain a rare event, which could not be 
considered a representative element of the source population of 
ultra-high energy cosmic rays. Overall, one finds that this source 
must contribute to a fraction 10 % of the flux at energies well be- 
low the GZK cut-off in order to be detectable in gamma-rays. It 
is important to note that this conclusion does not depend on the 
degree of beaming of the emission of the source, as we discuss 
in SectionlH 

To conclude this section on a more optimistic note, let us 
discuss an effect that could improve the prospects for detec- 
tion. In the above di s cussio n and in our simulations, following 
iGabici & AharonianI (l2005h . we have assumed that all of syn- 
chrotron emission takes place in the immediate environment of 
the source, as the magnetic field drops to low values outside the 
filament, where inverse Compton cascades would therefore be- 
come the dominant source of energy loss of the pairs. However, 
since the typical energy loss length of ultrahigh energy cosmic 
rays exceeds the typical distance between two filaments, one 
must in principle allow for the possibility that the beam of ultra- 
high energy protons seeds secondary pairs in intervening magne- 
tized filaments, which would then also contribute to the observed 
synchrotron emission. During the crossing of a filament, ultra- 



high energy protons are deflected by an angle 66i <s; 1 (unless 
the primary cosmic ray carries a large charge or the filament is of 
unusua l magnetisation, see the discussion in lKotera & Lemoiii3 
'2QP8b'). Then, as long as the small deflection limit applies, the 
problem is essentially one-dimensional and each filament con- 
tributes to the gamma ray flux a fraction comparable to that of 
the source environment. The overall synchrotron signal would 
then be amplified by a factor /l//f, where A denotes the relevant 
energy loss length, /I ~ 100 Mpc for pion production of 10^° eV 
protons, A ~ \ Gpc for pair production of 10''' eV protons and 
If ~ 40 Mpc is the typical separation between two filaments. 
Overall, one might expect an amplification as large as an order 
of magnitude for the contribution to the flux due to pair produc- 
tion, for sources located at several hundreds of Mpc, and a factor 
of a few for sources located at 100 Mpc; and, correspondingly, 
an amplification factor ~ 2-3 for the contribution to the flux due 
to pion production. Since the pair production channel contributes 
~ 10% of the total flux at fi^ ~ 10-100 GeV for a pure pro- 
ton composition, the amplification factor should remain limited 
to a factor of a few, unless the interaction length with magne- 
tized regions is substantially smaller than the above (expected) 
If =40 Mpc. 



3.3. Inverse Compton cascades 

Let us briefly discuss the gamma ray signal expected from 
Compton cascades of ultra-high energy photons and pairs in- 
jected in the intergalactic medium. The physics of these cas - 
cades has been discussed in detail in 'Wdowczvk et al.' (1972^ 
Prothe r oe (1986), Pro theroe & Stan^v (1993), Aharonian et^] 
(11994 . iFerrigno et al.l (l2004 and lGabici & AharonianI tOQlh . 
The angular extent and time distribution of GeV-TeV gamma- 
rays resulting from inverse Compton cascades seeded by 
ultra-high energy cosmic rays produced by gamma-ray bursts 
have also been discussed by WMiMn_& Coppi ( 1996). Inverse 
Compton cascades in t he steady state re g ime h ave been con- 
sidered in the study of lArmengaud et aP (l2006h (for a source 
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smaller than 10'^ eV, so that the contribution of lower energy 
particles can be neglected in a first approximation. For photo- 
pair production, the fraction transfered is Xe.ee - d/l Gpc of 
Lgjc) - La-{> 10'^ eV) up to d ~ 1 Gpc. For pion production, 
the fraction of energy transfered is roughly j^^^^^j ^ d/lOO Mpc of 
Lcr(> 6 10'^ eV) in the continuous energy loss approximation. 
At distances 100 Mpc < d < 1 Gpc, the fraction of L^^ 19 in- 
jected into secondary pairs and photons thus ranges from ~ 0.5 
for d = 100 Mpc to ~ 1 at = 1 Gpc; in short, it is expected 
to be of order unity or slightly less. All the energy injected in 
this way in sufficiently unmagnetized regions (see below) will 

be deposited through the inverse Compton cascade in the sub- 

1 y 

TeV range, with a typical energy flux dependence oc Ey" up 
to some maximal energy E 
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Fig. 8. Angular profiles of the images of the sources represented 
in Fig.|2l We represent the gamma ray flux integrated over ener- 
gies Ey = 1-100 GeV averaged over angular bins, for a filament 
seen along its axis, at 1 Gpc and Le,i9 - 10*^ erg s"' (black solid 
line), and at 100 Mpc and 19 = lO'*"* erg s"' (green dashed 
line). The black stars and green crosses present the correspond- 
ing integrated flux up to a given angular extension in the sky 9. 



located in a cluster of galaxies) but dismissed in the study of 
iGabici & AharonianI (l2005i) because of the dilution of the emit- 
ted flux through the large deflection of the pairs in the low energy 
range of the cascade. Indeed, the effective inverse Compton cool- 
ing length of electrons of energy Eg < 100 TeV can be written as 
Xey - 3.5 kpc (Zip/ 100 TeV)"' and on this distance scale, the de- 
flection imparted by a magnetic field of coherence length Ag » 

reads 6», ~ x,y/rL,e ~ 3 x lQ-\Ee/10QTeWr\B/lO-^-G). 
Then, assuming that the last pair of the cascade carries an en- 
ergy Ef\„ ~ 20 TeV (so that the photon produced through the 
interaction with the CMB carries a typical energy < 1 TeV), 
one finds that a magnetic field larger than ~ 10 G isotropizes 
the low energy cascade, i n agreement with the estimates of 
iGabici & AharonianI (l2005h . 

This situation is modified when one takes into account the 
inhomogeneous distribution of extra-galactic magnetic fields, as 
we now discuss. Primary cosmic rays, upon traveling through 
the voids of large scale structure may inject secondary pairs 
which undergo inverse Compton cascades in these unmagnetized 
regions. If the field in such regions is smaller than the above 
10"'^ G, then the cascade will transmit its energy in forward 
<TeV photons. Of course, depending on the exact value of B 
where the cascade ends, the resulting image will be spread by 
some finite angle. Since we are interested in sharply peaked im- 
ages, let us consider a typical angular size 6 and ignore those 
regions in which the magnetic field is large enough to give a 
contribution to the image on a size larger than 9. For 9 I, 
the problem remains one-dimensional as before, and one can 
compute the total energy injected in inverse Compton cascades 
within 9, as follows. 

The luminosity injected in secondary pairs and photons up 
to distance d is written XeLa-(> E). Since we are interested 
in the signatures of ultrahigh cosmic ray sources, we require 
that E > 10'"^ eV; for protons, the energy loss length due 
pair production moreover increases dramatically as E becomes 



y.max 



1 - 10 TeV beyond which 
the Universe is opaque to gamma rays on the distance scale d 
(Ferrigno et al. 2004). Neglecting any redshift dependence for 
simplicity, the gamma-ray energy flux per unit energy interval 
may then be approximated as: 



'''ae: 



/id(< Be)Xe 



1/2 



87r£/2 \Ey 

^ 2.5 X 10""' GeV cm"^ s"' /id(< Bg)xe 
Lea9 I d \ ^ I E^ 



1042 erg/s \ 100 Mpc 



1/2 



(6) 



where /id(< Be) denotes the one-dimensional filling factor, i.e. 
the fraction of the line of sight in which the magnetic field is 
smaller than the value Bg such that the deflection of the low en- 
ergy cascade is 9. For reference. Be ~ 2 x lO^'^G for 9 ~ 1°. 
In general, one finds in the literature the three-dimensional fill- 
ing factor /3d, but /id(< Bg) ~ /3d(< Bg) up to a numerical 
prefactor of order unity that depends on the geometry of the 
structures. Interestingly enough, the amount of magnetization of 
the voids of large scale structure is directly related to the origin 
of large scale magnetic fields. Obviously, if galactic and clus- 
ter magnetic fields originate from a seed field produced in a 
homogeneous way with a present day strength B » lO '^G^ 
then the above gamma ray flux will be diluted to large angular 
scales, hence below detection threshold. However, if the seed 
field, extrapolated to present day values is much lower than 
this value, or if most of the magnetic enrichment of the in- 
tergalactic medium results from the pollution by star forming 
galaxies and radio-galaxies, then on e should expec t /id(< Bg) 
to be non negligible. For instance, iDonnert et al.1 (I2009h ob- 
tain /3d(< lO '^G) ~ 0(1) in such models. Given the sensi- 
tivity of current gamma ray experiments such as H.E.S.S. in 
the TeV energy range, the inverse Compton cascades might 
then produce degree-size detectable halos for source luminosi- 
ties > 1043[/3d(< 10-14 G)]-i(^^/i00Mpc)-2erg/s. For future in- 
struments such as CTA, the detectability condition will be of 
order Ls > 10'''[/3d(< 10"'4Q)]-i(^/ioOMpc)--erg/s. The ex- 
pected flux level remains quite uncertain as it depends mostly on 
the configuration of the extragalactic magnetic field in the voids, 
contrarily to the synchrotron signal which is mainly controlled 
by the source luminosity. 

We note that intergalactic magnetic fields of strength B < 
10 G might be probed through the delay time of the high en 



ergy afterglow of gamma-ray bursts (|Plagal[T995 : Ichiki et al 
2008h or flie GeV e mission around blazar s jAharonian et al 



1994- 'Dai et al."2002'; 'd'A vezac et all 120071: lElviv et al.ll2009[ 



Nerongy & Semikoz 20Q9), as discussed earlier. As the present 
paper was being refereed, some first estimates on lower bounds 
for the average magnetic field appeared, giving B > 10 - 
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10" '^ G. see | Neronov & Vov^ (l2010l) . lAndo & Kusenkol (l2010l) 
and lDolag et al.l(l2010l) . 



3.4. Cen A 

In view of the results of section 13.11 one may want to consider 
also the case of mildly powerful but nearby sources. One such 
potential source is Cen A, which has attracted a considerable 
amount of attention, as it is the nearest radio-galaxy (3.8 Mpc) 
and more recently, because a fraction of the Pierre Auger events 
above 60 EeV have c lustered within 10 - 20° of this source. As 
discussed in detail in lLemoine & WaxmanI (l2009l) . this source is 
likely too weak to accelerate protons to > lO'^ eV in steady state, 
but heavy nuclei might possibly be accelerated up to GZK ener- 
gies. The acceleration of protons to ultrahigh energies in flaring 
episod es of high luminosity has been proposed in Dermer et al,! 
(l2009l) . It is also important to recall that the apparent cluster- 
ing in this direction can be naturally explained thanks to the 
large concentration of matter in the local Universe, in the di- 
rection of Cen A but located further away. As discussed in 
iLemoine & Waxrnan ( 2009), some of the events detected toward 
Cen A could also be attributed to rare and powerful bursting 
sources located in the Cen A host galaxy, such as gamma-ray 
bursts. Due to the scattering of the particles on the magnetized 
lobes, one would not be able to distinguish such bursting sources 
from a source located in the core of Cen A. 

We have argued throughout this paper that a clear signature 
of ultrahigh energy cosmic ray propagation could be produced 
in presence of an extended and strong magnetized region around 
the source. The magnetized lobes of Cen A provide such a site, 
with an extension of Riobe ~ 100 kpc, and a magnetic field of av- 
erage intensity B]n Uf ~ 1 /"G and coherence length /liobe ~ 20 kpc 
dFeain et al.ll2009h . In terms of gamma ray fluxes, one can calcu- 
late however that the only advantage Cen A presents compared 
to average types of sources described in section 13.11 lies in its 
proximity. As noted before, the bulk of gamma ray emission will 
be due to electrons and positrons produced by pion production 
with Ef > 10'^ eV. The flux of these secondary electron/positron 
emission scales as the distance r the primary protons travels 
through, as long as r <k /I, with A the energy loss length by 
pion production. In our situation, in order to calculate the syn- 
chrotron gamma ray flux, this distance corresponds to the size 
of the magnetized structure with B > 0.5 nG(£'J10''^ eV)"^^/"*, 
namely r ~ Riobe in our toy-model. One may also note that con- 
sidering the magnetic field strength in the lobes, the gamma ray 
flux will peak around fiysyn ~ 6.8 TeV Bio\,gE^ jg, according to 
Eq.[2] Scaling the expected flux from the lobes around 10 TeV 
^lobe.iOTeV to the flux that we calculated in section [TTI for stan- 
dard sources embedded in filaments around 10 GeV /^fii,ioGeV, 
one obtains: 



-Flobe.lOTeV 



da 
2x 10 



Ru 



lO^^erg/s 5 Mpc 
GeV s^'cm"" , 



fil.lOGeV 



(7) 
(8) 



which is far below current instrument sensitivities. Here, dcenA - 
3.8 Mpc and da = 100 Mpc are the distance of Cen A and of 
our model filament respectively, and LcenA is the cosmic ray lu- 
minosity of Cen A which we set to 10^^ erg s"' according to 
constraints f rom ultrahigh energy co smic ray observations (see 
for example iHardcastle et al.ll2009h . Moreover, Cen A covers 
^source ~ 10° in the sky, inducing a sensitivity loss of a factor 
^source /^^PSF ~ 10 for the Fermi telescope (ffpsp represents the 
angular resolution of the instrument). 



All in all, taking into account all these disadvantages, we 
conclude that the synchrotron emission produced in the lobes 
of Cen A would be three orders of magnitudes weaker than the 
flux obtained for the average types of sources described in sec- 
tion O 

Our numerical estimate using the method introduced in sec- 
tion |2]2] and modeling the lobes as spherical uniformly magne- 
tized structures confirms this calculation. We note that the con- 
tribution of the pair production channel is not as suppressed 
as the other components. Due to the strength of the magnetic 
field, the electrons and positrons radiating in synchrotron around 
0.1-1 GeV are of lower energy than for the standard filament 
case {Ef ~ 10 MeV) and thus are more numerous. At these ener- 
gies, the synchrotron emission in the filament was also partially 
suppressed by the shortening of the inverse Compton loss length 
as compared to the synchrotron loss length. Because of the high 
magnetic field intensity, this effect is less pronounced in the lobe 
case. 

Recently Fermi LAT discovered a gam ma ray emission ema - 
nating from the giant radio lobes of Cen A ( ICheung et al.ll2010h . 
In light of our discussion, the ultrahigh energy cosmic ray ori- 
gin of this emission can be excluded. It could rather be inter- 
preted as an inverse Compton scattered relic radiation on the 
C MB and infrared and optical photon backgrounds, as proposed 
bv lCheungetaP (l2010l) . ^ 

Let us note additionally that iKacheMeB et alj ( l2009h have 
provided estimates for gamma ray fluxes emitted as a result 
of the production of ultrahigh energy cosmic rays in Cen A. 
These authors conclude that a substantial signal could be pro- 
duced, well within the detection capabilities of Fermi and CTA. 
However this signal results from interactions of cosmic rays ac- 
celerated inside the source, not outside. As we discussed before, 
such a gamma ray signal, i.e. that could not be angularly re- 
solved from the source, could not be interpreted without ambi- 
guity. In particular, one would not be able to discriminate it from 
a secondary signal produced by multi-TeV leptons or hadrons 
accelerated in the same source. Furtherm ore, the gamma ray 
signal discussed in IKacheMeB et all ( l2009h is mainly attributed 
to cosmic rays at energies < 10'^ eV, not to ultrahigh cosmic 
rays above the ankle, although the generation spectrum in the 
source is normalized to the Pierre Auger Observatory events at 
> 6x10''^ eV. The conclusions thus depend rather strongly on the 
shape of the injection spectrum as well as on the phenomenolog- 
ical modelling of acceleration, as made clear by Fig. 1 of their 
study. 

Finally one may turn to the highest energies and search for 
ultrahigh energy cosmic ray signatures from Cen A by looking 
directly at the ultrahigh energy secondary photons produced dur- 
ing the propagation. These phot ons should indeed n ot experience 
cascading at this close distance. iTavlor et"al1 ( |2009|) have studied 
this potential signature and concluded that Auger should be able 
to detect 0.05 - 0.075 photon per year from Cen A, assuming 
that it is responsible for 10% of the > 6 x 10''' eV cosmic ray 
flux, and assuming a 25% efficiency for photon discrimination. 

4. Discussion and conclusions 

Let us first recap our main findings. We have discussed the de- 
tectability of gamma ray halos of ultrahigh energy cosmic ray 
sources, relaxing most of the assumptions made in previous stud- 
ies. We have focused our study on the synchrotron signal emitted 
by secondary pairs, which offers a possibility of unambiguous 
detection. In this study, we have taken into account the inhomo- 
geneous distribution of the magnetic field in the source environ- 
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ment and examined the effects of non pure proton composition 
on the gamma ray signal. We find that the gamma ray flux predic- 
tions are rather robust with respect to the uncertainties attached 
to these physical parameters. The normalization and hence the 
detectability of this flux ultimately depends on the energy in- 
jected in the primary cosmic rays, for realistic values of the mag- 
netic field. We have further demonstrated that the average type of 
sources contributing to the ultrahigh energy cosmic ray spectrum 
produces a gamma ray flux more than two orders of magnitudes 
lower than the sensitivity of the current and upcoming instru- 
ments. The case of rare powerful sources contributing to more 
than 10% of the cosmic ray flux is far more promising in terms 
of detectability. We have found that gamma ray signatures of 
those sources could be detectable provided that they are located 
far enough not to overshoot the observed cosmic ray spectrum. If 
the extended emission of such signatures were resolved (which 
should be the case with Fermi and CTA), such a detection would 
provide a distinctive proof of acceleration of cosmic rays to ener- 
gies > lO'^ eV. However, as we have discussed, such gamma ray 
sources should remain exceptional and could not be considered 
as representatives of the population of ultrahigh energy cosmic 
ray sources. We have also discussed the deflection and the di- 
lution of the Compton cascading gamma ray signal at sub-TeV 
energies. In this case, the flux prediction depends sensitively on 
the filling factor of the magnetic field with values smaller smaller 
than ~ 10"'"* G in the voids of large scale structure. If this latter 
is larger than a few percents, as suggested by some scenarios of 
the origin of cosmic magnetic fields, then the inverse Compton 
cascade could provide a degree size image of the ultra-high cos- 
mic ray source, with a flux slightly larger than that associated to 
the synchrotron component. 

Finally, we have also discussed the detection of nearby 
sources, considering the radiogalaxy Centaurus A as a protoypi- 
cal example. We have demonstrated that, //"ultrahigh energy cos- 
mic rays are produced in Cen A, their secondary gamma ray halo 
signature should not be detectable by current experiments. Our 
conclusion indicates in particular that the gamma ray emission 
detected by Fermi L AT from the giant radio lobes of Cen A 
(ICheung et al.ll2010l) cannot result from the synchrotron radia- 
tion of ultrahigh energy secondaries in that source. 

We need to stress that, in our discussion, we have system- 
atically assumed that the source was emitting cosmic rays at a 
steady rate with isotropic output. Let us discuss these assump- 
tions. If the source is rather of the flaring type, then, at a same 
cosmic ray source luminosity, the flux predictions will be re- 
duced by the ratio of the source flare duration Af to the dispersion 
in arrival times cr, that results from propagation in intergalactic 
magnetic fields. In theoretical models, Ar can take values rang- 
ing from a few seconds for gamma -ray bursts dMilgrom & Usovl 
Il996t IWaxmanl[T99l IVietrilll995h to days in the case of blazar 
flares, while cr, is typically expected to be o f the order of 10^ - 
10^ yrs for protons of energy > 70 Be V (see iKotera & Lemoin3 
l2008bl for a detailed discussion). In both cases, the large lumi- 
nosity during the flaring episodes is not sufficient to compen- 
sate for the dilution in arrival times. Of course, if the source is a 
repeater, with a timescale between two flaring episodes shorter 
than the intergalactic dispersion timescale, then one recovers a 
steady emitting situation, with an effective luminosity reduced 
by the fraction of time during which the source is active. 

This latter situation may well apply to powerful radio- 
galaxies that produce ultra-high energy cosmic rays with flaring 
activity, since the size and magnetization of the jets and lobes are 
sufficient to deflect ultra-high energy cosmic rays by an angle of 
order unity, theredy dispersing them in time over ~ R\ohe/c ~ 



3x10^ yr, with /?iobe ~ 100 kpc the typical transverse scale of the 
lobes. Even if ultra-high energy cosmic rays are ejected as neu- 
trons, as discussed for instance by Dermer et al. (2009), a sub- 
stantial fraction of these cosmic rays will suffer dispersion as the 
decay length of neutrons of energy E„ is 0.9(£'„/10^°eV)Mpc. 
In such a configuration, the source would effectively emit its cos- 
mic rays isotropically, as we have assumed. 

Let us nonetheless discuss the case of beamed sources, for 
the sake of completeness. If one considers a single source, emit- 
ting cosmic rays in a solid angle AOs < 47r, then at a same to- 
tal luminosity as an isotropic source, the luminosity emitted per 
steradian is larger by 47r/AOs and so is the expected gamma-ray 
flux (in the limit of small deflection around the source, which 
is necessary for detection). In other words, the required total 
cosmic-ray luminosity per source for detection of gamma rays is 
lowered by a factor Afis/(47r). However, at a fixed contribution 
of the source to the cosmic ray flux, the expected gamma-ray flux 
does not depend on the beaming angle. This is because both the 
cosmic ray flux and the gamma-ray flux are determined by the 
number of visible sources times the luminosity per steradiarQ In 
particular, our conclusion that a source must contribute to more 
than 10% of the ultra-high energy cosmic ray flux in order for 
the synchrotron signal to be detectable remains unchanged if 

AOs < 47T. 
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